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Abstract: The dynamical responses of ferromagnet to the propagating elec- 
tromagnetic field wave passing through it are modelled and studied here by 
Monte Carlo simulation in two dimensional Ising ferromagnet. Here, the 
electromagnetic wave is linearly polarised in such a way that the direction 
of magnetic field is parallel to that of the magnetic momemts (spins). The 
coherent propagating mode of spin-clusters is observed. The time average 
magnetisation over the full cycle (time) of the field defines the order param- 
eter of the dynamic transition. Depending on the value of the temperature 
and the amplitude of the propagating magnetic field wave, a dynamical phase 
transition is observed. The dynamic transition was detected by studying the 
temperature dependences of the dynamic order parameter, the variance of the 
dynamic order parameter, the derivative of the dynamic order parameter and 
the dynamic specific heat. The phase boundaries of the dynamic transitions 
were drawn for two different values of the wave lengths of the propagating 
magnetic field wave. The phase boundary was observed to shrink (inward) for 
lower speed of propagation of the EM wave. The divergence of the releavant 
length scale was observed at the transition point. 
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I. Introduction. 

The dyamical response of Ising ferromagnet to a time dependent magnetic 
field has become an active field of research [1]. The hysteretic responses and 
the nonequilibrium dynamic phase transitions are two major foci of atten- 
tion. The scaling behaviour [2] of hysteretic loop area with the amplitude, 
frequency of the sinusoidally oscillating magnetic field is the main outcome of 
the research. Another interesting aspect is the nonequilibrium dynamic phase 
transition which yields variety of interesting nonequilibrium behaviours and 
promted the researchers to take continuous attention in this field. Histor- 
ically, some important observations like (i) divergences of dynamic specific 
heat and relaxation time near the transition point [3], (ii) divergence of the 
scale of length near the transition point [4], (iii) studies regarding the exis- 
tence of tricritical point [5, 6], (iv) the relation with the stochastic resonance 
[5] and the hysteretic loss [7], enriched the field and established the transition 
as a thermodynamic phase transition. Very recently, a surface dynamic phase 
transition [8] is observed in kinetic Ising ferromagnet driven by oscillating 
magnetic field. The dynamic phase transition was studied experiment ally [9] 
also, in ultrathin Co film on Cu(100) system by surface magneto-optic Kerr 
effect. 

This dynamic phase transition is also observed in other magnetic mod- 
els. The dynamic phase transition was observed [10] in anisotropic clas- 
sical Heisenberg model and in XY model [11]. The multiple (surface and 
bulk) dynamic transition was observed [12] in Heisenberg model (with bi- 
linear exchange anisotropy) and it was found [13] also in Heisenberg ferro- 
magnet driven by polarised magnetic field. The dynamic transition were 
observed [14] in kinetic spin-3/2 Blume-Capel model and in Blume-Emery- 
Griffith model[15] studied by meanfield approximation. The dynamic phase 
transition was studied by Monte Carlo simulation [16] and by meanfield 
calculation[17] in Ising metamagnet. 

However, all the studied mentioned so far, were done by sinusoidally oscil- 
lating magnetic field which was uniform over the space (lattice) at any instant 
of time. In those studies, the spatio-temporal variation of applied magnetic 
field was not considered. One such spatio-temporal variation of applied mag- 
netic field would be the propagating magnetic field wave. In reality, if the 
electromagnetic wave passes through the ferromagnet, the spatio-temporally 
varying magnetic field coupled with the spin, would affect the dynamic na- 
ture of the system. Here also dynamic transition would be observed. Very 
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recently, it is reported briefly [18] that propagating magnetic field wave would 
lead to dynamical phase transition in Ising ferro magnet. A pinned phase and 
a phase of coherent motion of spin-clusters were observed recently [19] in ran- 
dom field Ising model swept by propagating magnetic field wave. Here the 
nonequilibrium dynamic phase transition is athermal and tunned by random 
disorder. A rich dynamical phase boundary was also drawn. 

Here in this paper, the nonequilibrium dynamic phase transition is stud- 
ied extensively in two dimensional Ising ferromagnet swept by propagating 
electromagnetic field wave. The technique employed here is Monte Carlo 
simulation. The phase boundary of the dynamical transition is drawn in this 
study. The paper is organised as follows: The model and the simulation 
technique are discussed in section-II, the numerical results are reported in 
section-Ill and the papers ends with a summary, in section-IV. 

II. Model and Simulation. 

The Hamiltonian (time dependent) representing the two dimensional Ising 
ferromagnet in presence of a propagating electromagnetic field wave (having 
spatio-temporal variation) can be written as 

H{t) = -JEs(x, y, t)s(x', y', t) - T,h(x, y, t)s(x, y, t) (1) 

The s(x,y,t) represents the Ising spin variable (±1) at lattice site (x,y) at 
time t on a square lattice of linear size L. J(> 0) is the ferromagnetic (taken 
here as uniform) interaction strength. The first sum for Ising spin-spin inter- 
action is carried over the nearest neighbours only. The h(x, y, t) is the value 
of the magnetic field (at point (x, y) and at any time t) of the propagating 
electromagnetic wave. It may be noted here that the electroagnetic wave is 
linearly polarised in such a way that the direction of magnetic field is parallel 
to that of the spins. For a propagating magnetic field wave h(x, y, t) takes 
the form 

h(x, y, t) = h cos(2irft - 2ny/\) (2) 

The ho, f and A represent the amplitude, frequency and the wavelength 
respectively of the propagating electromagnetic field wave which propagates 
along the y-direction. In the present simulation, alxL square lattice (in 
contact with a heat reservoir of temperature T) is considered. The boundary 
condition, used here, is periodic in both (x and y) direction. The initial (t = 
0) configuration, as half of the total number (selected randomly) of spins are 
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up (s(x,y,t = 0) = +1 for all x and y)), is taken here. This configuration of 
spins, corresponds to the high temperature disordered (paramagnetic) phase. 
The spins are updated randomly (a site (x, y) is chosen at random) and the 
spin flip occurs (at temperature T) according to the Metropolis rate [20] W 

W(s -)■ -s) = Min [exp (—AE/kT), 1], (3) 

where AE is the change in energy due to spin flip and k is Boltzmann con- 
stant. L 2 such random updates of spins constitutes the unit time step here 
and is called Monte Carlo Step (MCS). Here, the value of magnetic field is 
measured in the unit of J and the temperature is measured in the unit of 
J/k. The dynamical steady state , at any temperature, is reached by cooling 
the system slowly (in small step of temperature, AT = 0.02) from the high 
temperature configuration. The frequency of the propagating magnetic field 
was taken / = 0.01 throughout the study. The total length of the simulation 
is 2 x 10 5 MCS and first 10 5 MCS transient data ware discarded. The data 
are taken by averaging over 10 5 MCS. Since the frequency of the propagating 
field is f — 0.01, the complete cycle of the field requires 100 MCS. So, in 10 5 
MCS, 10 3 numbers of cycles of the propagating EM wave are present. The 
time averaged data over the full cycle (100 MCS) of the propagating field are 
further averaged over 1000 cycles. 

III. Results. 

In this study, a square lattice of size L = 100 is considered. The steady 
state dynamical behaviours of the spins are studied here. The amplitude, 
frequency and the wavelength of the propagating wave are taken ho = 0.6, 
/ = 0.01 and A = 25 respectively. The magnetic field is propagating along the 
y-direction (vertically upward in the graphs shown here). The temperature 
of the system is taken T = 1.50. The configuration of the spins at any 
instant of time t = 100100 MCS is shown in Fig-l(a). Here, it is noted 
that, the clusters of spins (down) are formed in strips and these strips move 
coherently as time goes on. The propagation of the spin-strips are clear in 
Fig-l(b), where the snapshot was taken at later instant t = 100125. The 
similar study was performed at a lower temperature T = 1.26 (with all other 
parameters of the propagating field remain same). Here, the spin clusters are 
observed to form but not in the shape of strips (as observed in the case of 
higher temperature T = 1.50 , mentioned above). This is shown in Fig- 1(c), 
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at any instant t = 100100 MCS. These irregularly shaped spin-clusters are 
observed to propagate (along the direction of propagating magnetic field), 
which is clear from Fig-l(d) (for t = 100125 MCS). 

To show the propagations of these spin-clusters, the instantaneous line 
magnetisation m(y,t) = J s(x,y,t)dx/L was plotted against y at any par- 
ticular instant t = 100100 MCS. This is shown in Fig-2(a) (compare with 
Fig-l(a)). The periodic variation of m(y,t) along y-direction is found. This 
was observed to propagate (see Fig-2(b) and compare with Fig- 1(b) when 
shown at later time t = 100125 MCS. It may be noted here, that the line 
magnetisation is periodic (with y) at any instant of time t. This is also peri- 
odic in time t at any position y. In both the cases, the oscillation is symmetric 
about zero (for higher temperature T = 1.50). Here, the time average mag- 
netisation over a full cycle of the propagating field is Q = J J m(y, t)dydt, 
becomes zero (due to symmetric oscillation about zero). This is dynamically 
symmetric disordered (Q = 0) phase. 

Now, for lower temperature T = 1.26, the spatio-temporal periodicity is 
lost. The symmetric-oscillation (about zero) is lost here. As a consequence, 
the time average magnetisation over a full cycle of the propagating field, be- 
comes nonzero (Q ^ 0). These are shown in Fig-2(c) and Fig-2(d) (compare 
with Fig-l(c) and Fig-l(d) respectively). But the spin-clusters propagates, 
in this case. So, as the temperature decreases, Q becomes nonzero (lower 
temperature) from a zero value (higher temperature). This is dynamically 
symmetry-broken disordered (Q ^ 0) phase. 

The temperature variations of the dynamic order parameter Q, its stan- 
dard deviation < (SQ) 2 > are studied. The dynamic energy is E = ^- § H (t)dt 
and the dynamic specific heat is C = 4p. The derivatives are calculated nu- 
merically by using the three points central difference formula[21] All these 
quantities are calculated statistically (averaged) over 1000 different values. 
The temperature variations of Q, < (5Q) 2 > and C are studied fow two 
different values of the amplitude of the propagating magnetic field and are 
shown in Fig-3. As the temperature decreases Q starts to grow from zero 
and near the transition point T d (depending upon the value of ho and A) it 
becomes nonzero. Near the transition temperatures the < (5Q) 2 > and C 
show shap peaks and show a sharp dips. From the figure it is also evi- 
dent that the transition occurs at lower temperature (T^) for higher values 
of the field amplitude (ho)- In this case, for A = 25, the transitions occur at 
T d = 1.90 and T d = 1.28 for h = 0.3 and h = 0.6 respectively. These values 
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of the transition temperatures are obtained from the position of sharp dip of 
the and corresponding sharp peaks of < (5Q) 2 > and C shown in Fig-3. 
Collecting all values of the transition temperatures (depending on the values 
of h ), the dynamical phase boundary is obtained. 

This dynamic transition temperature (7^) was observed to depend on the 
wave length (A) of the propagating magnetic field wave. The temperature 
dependences of Q, < (5Q) 2 > and C are studied and shown in Fig-4, 
for two different values of A (= 25 and 50). From the figure, it is clear that, 
transition occurs at higher temperature for higher value of the wavelength 
(A). To be precise, for ho = 0.3, the transitions occur at Td = 1.88 and 
Td = 1.94 for A = 25 and A = 50 respectively. Here also, the values of the 
transition temperatures are obtained from the position of sharp dip of the 

and corresponding sharp peaks of < (5Q) 2 > and C (shown in Fig-4). So, 
the dynamical phase boundary should shift depending on the value of A. 

In Fig-5, the dynamical phase boundaries are drawn for two different 
values of A (=25 and 50), in the plane formed by Td and ho- It is observed that 
the boundary shrinks inward as the wavelength of the propagating magnetic 
field decreases. Since, the speed of propagating EM wave is c = /A, for 
fixed /(= 0.01 here), the phase boundary essentially shrinks inward (low 
teperature) for lower speed of propagation. 

The dynamic phase transition mentioned above is associated with the 
divergences of releavant length scale. For this reason the L 2 < (5Q 2 ) > 
is studied as a function of temperature T. Is is found that the peak of 
L 2 < (5Q) 2 > (observed at T d ) increases as L increases. This is shown in 
Fig-6. This result is quite conclusive to say that there exists the diverging 
length scale associated with the dynamic phase transition. It may be noted 
here, that this method was successfully employed [4] to show the diverging 
length scale associated to the dynamic transition in Ising ferromagnet driven 
by oscillating (but not propagating) magnetic field. 

IV. Summary. 

The dynamical responses of a ferromagnet to a polarised electromagnetic 
wave is modelled and studied here by Monte Carlo simulation in two di- 
mensional Ising ferromagnet. In the steady state, the coherent motion (in 
propagating mode) of spin clusters was observed in dynamically symmetric 
disordered phase. The time average magnetisation over the full cycle of the 
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propagating EM wave is a measure of a order parameter in the dynamic 
phase transition observed here. The dynamic phase transition observed here 
is of continuous type and found to be dependent on the amplitude and the 
wave length of the propagating polarised EM wave. Hence, a phase boundary 
(transition temperature as a function of the amplitude) is drawn for two dif- 
ferent values of the wavelengths of EM wave. The phase boundary was found 
to shrink (towards the lower temperature) for lower speed of propagation of 
EM wave. 

Another important observation is the divergence of releavant length scale 
near the dynamic transition. This observation supports the dynamic transi- 
tion to become a member of the well known equilibrium critical phenomenon. 
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Fig-1. The motion of spin-clusters of down spins (shown by black dots), 
swept by propagating magnetic field wave, for different values of (a) Time 
= 100100 MCS, T=1.5 and h = 0.6 (b) Time = 100125 MCS, T = 1.5 and 
h = 0.6 (c) Time = 100100 MCS, T = 1.26 and h = 0.6 (d) Time = 100125 
MCS, T = 1.26 and h = 0.6. 
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Fig-2. The propagation of field (x) and the line magnetisation (+) for 
various values (from top to bottom) of (a) Time = 100100 MCS, T=1.5 
and ho = 0.6 (b) Time = 100125 MCS, T = 1.5 and h = 0.6 (c) Time = 
100100 MCS, T = 1.26 and h = 0.6 (d) Time = 100125 MCS, T = 1.26 and 
hn = 0.6. 



11 




dT -1U 

-20 



T 



Hlii 



x< 



(b) 



0.04 



< {5Q? % x 



T 



1 1 X 

X 


1 1 1 

, 1 (c) = 







c 



T 




Fig-3. The temperature (T) dependence of (a) Q, (b) jg, (c) < (5Q) 2 > 
and (d) C, for two different values of ho for propagating magnetic field wave 
having / = 0.01 and A = 25. In each figure, h = 0.3(+) and h = 0.6(x). 
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Fig-4. The temperature (T) dependence of (a) Q, (b) jg, (c) < (5Q) 2 > 
and (d) C, for two different values of A for propagating magnetic field wave 
having / = 0.01 and h = 0.3. In each figure, A = 25(+) and A = 50(x). 
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Fig-5. The phase diagram for dynamic phase transition by propagating 
magnetic field wave for two different values of wavelengths, A = 25(x) and 
A = 50(+). Here, / = 0.01. 
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Fig-6. The plot of temperature (T) versus L 2 < (5Q) 2 > for different system 
sizes. Here, h = 0.6, A = 25 and / = O.Of . 
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